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COUPLINGS OF BROWNIAN MOTIONS OF DETERMINISTIC DISTANCE IN MODEL 

SPACES OF CONSTANT CURVATURE 

MIHAI N. PASCU AND lONEL POPESCU 


Abstract. We consider the model space of constant curvature K and dimension n > 1 (Euclidean 
space for A = 0, sphere for K > 0 and hyperbolic space for K < 0), and we show that given a function 
p : [0, oo) —>■ [0, oo), p(0) = d{x, y) there exists a co-adapted coupling {X{t), Y(t)) of Brownian motions on 
such that p{t) — d{X (t), Y(t)) for every t > 0 if and only if p is continuous and satisfies for almost every 
i > 0 the differential inequality 

-(n - l)^tan (^) < p'(t) < -(n - 1)7^tan (^) + 

In other words, we characterize all co-adapted couplings of Brownian motions on the model space for 
which the distance between the processes is deterministic. In addition, the construction of the coupling is 
explicit for every choice of p satisfying the above hypotheses. 


1. Introduction 

The initial motivation for writing the present article was to investigate a stochastic version of the 
celebrated "Lion and Man" problem of Rado ([17]), which asks for winning strategies in a game in which 
a Man is chased by a Lion in a circular arena (to make the game interesting, it is assumed that both 
trajectories are smooth and of unit speed). Substituting the circular arena with a model space, namely 
a manifold of constant curvature, and the smooth trajectories of the Man and the Lion by Brownian 
trajectories on this space, the problem becomes a question about the existence of Brownian couplings for 
which the distance between the two processes is bounded away from zero (of particular interest being 
the case when it is constant), respectively when it approaches (or becomes) zero. 

In this paper the model spaces are classified by the curvature constant K. In the case of A = 0 this is 
the Euclidean space, in the case K > 0 this is the sphere of radius 1 / \/A while in the case of iT < 0 is the 
hyperbolic space with a certain metric. 

Regarding couplings, there are at least two cases of interest. The first is the case when one wishes 
to have a fast coupling time (the Lion's strategy), and the second is the case when one desires a slow 
(infinite) coupling time (the Man's strategy; think also of the case of a Brownian target being chased 
by a Brownian missile). There are various notions of fast/slow couplings in the literature: see, e.g., 
[8] for "optimal couplings", [7] for "efficient coupling", or [13] for "maximal couplings". In the case 
of Euclidean Brownian motions, the extremes of the coupling time are achieved by the mirror coupling, 
respectively by the synchronous coupling. 

The mirror coupling was first introduced in the Euclidean case by Lindvall and Rogers [16], and then 
extended to processes defined on manifolds by Cranston [9] and Kendall [14], the so-called Cranston- 
Kendall mirror coupling. For a recent extension of the mirror coupling to the case when the two (reflecting) 
Brownian motions live in different domains, see [19]. It turns out that this coupling is a very useful 
and versatile construction when it comes to various geometric and analytic properties on manifolds. 
For instance, it was shown in [14], for the case of manifolds with Ricci curvature bounded uniformly 
from below by a positive constant, that mirror coupled Brownian motions always meet in finite time. 
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Informally, the mirror coupling makes the Brownian motions move toward each other in the geodesic 
direction, and parallel to each other in the orthogonal direction. 

Under synchronous coupling, the Brownian motions move parallel to each other both in the geodesic 
direction and in the orthogonal direction. This coupling was used for example in [3] to prove the validity 
of Hot Spots conjecture for obtuse triangles, or in [2] to prove Harnack inequalities and heat kernel esti¬ 
mates on manifolds. On a different note, parametrized families of Brownian motions were constructed 
in [1] and [22]. 

In a different direction, the notion of shy coupling of Brownian motions was introduced in [4] and was 
further studied in [5] and [15]. It is a refinement of a slow coupling, for which with positive probability, 
the distance between the two processes stays positive for all times. A stronger version of shyness (e- 
shyness, e > 0) asserts that with positive probability, the distance between the processes is greater than 
e. In this paper we use this latter version of shyness, in the stronger sense that the distance between the 
processes is greater than e with probability 1. 

Thus, there are at least two cases of interest regarding Brownian couplings: on one hand is the case of 
couplings under which the Brownian motions get close to each other fast, and, on the other, is the case 
of couplings under which the particles eventually stay away from each other. 

Our focus in this paper is the following problem. Let p be a continuous function of time with p(0) being 
the initial distance between the Man and the Lion. Does the Brownian Man (or Lion) have a strategy in 
this game, such that the distance between him and the Lion (or Man) at time t is precisely p{t)? Our 
goal is to determine the class of functions p for which this is possible together with a description of the 
coupling. 

In the Euclidean case, given a function p : [0, oo) —> [0, oo) we show that there exists a co-adapted 
coupling {X{t),Y(t)) of Brownian motions such that d{X{t),Y{t)) = p{t) for all f > 0 if and only if p is 
continuous and satisfies a.e. the inequality 

0</(t)<!^. «>0. 

In particular the only co-adapted coupling with deterministic non-increasing distance function is the 
translation coupling. This is the content of Theorem 3.1. As a consequence, the extreme growth of the 
distance function is achieved for the distance function p{t) = \J p^(0) -|- 2(re — l)t. In terms of the stochas¬ 
tic version of the Lion and Man problem, this shows that there is no winning strategy for the Brownian 
Lion, while any co-adapted coupling with deterministic distance function is a winning strategy for the 
Brownian Man - the best choice is for the distance function which grows with square root in time. 

The other case study is that of the unit sphere C (n > 1), and the analogous result is presented 
in Theorem 4.1. The content is that there exists a co-adapted coupling of Brownian motions {X{t),Y(t)) 
on such that d{X{t), Y(t)) = p{t) for all t > 0 if and only if p : [0, cx)) —>• [0, cx)) is continuous and 
satisfies a.e. the differential inequality 

- 1) tan < p'(t) < - 1) tan t > 0. 

An important instance of this result is the case when p is a constant function, which gives the exis¬ 
tence of a fixed-distance coupling of spherical Brownian motions. The two other interesting particu¬ 
lar cases correspond to the cases of equality in the above inequality, being explicitly given by p{t) = 
2 arcsin sin(p(0)/2)), and p{t) = 2arccos (e“("“^)*/^cos(p(0)/2)). Both couplings are partic¬ 

ular cases of shy couplings: in the former case the two Brownian motions approach each other expo¬ 
nentially fast, but do not couple, and in the latter case the two Brownian motions repel each other at an 
exponential rate to the maximum distance allowed on the sphere (the processes become antipodal in the 
limit). In terms of the stochastic version of the Lion and Man problem, this shows that the Brownian 
Lion has a strategy which brings him within e > 0 from the Brownian Man, exponentially fast in time, 
and the Brownian Man has a strategy of safety, which increases his distance from the Brownian Lion 
to the maximum distance allowable on the sphere. Interestingly, the Brownian Man also has a strategy 
which keeps the Brownian Lion at fixed distance for all times, particularly frustrating from the point of 
view of the Brownian Lion. 
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The remaining case of a constant curvature manifold is that of the n-dimensional hyperbolic space H” 
(n > 1). In this case we show that there exists a coupling {X (t), Y(t)) of Brownian motions on H" with 
d{X{t),Y(t)) = p{t) if and only if p : [0, oo) —)• [0, oo) is a continuous function which satisfies almost 
ever}rwhere the inequality 

(« - Dtanh < p'(t) < („- l)tanh t > 0. 

From the point of view of the Brownian Lion and the Brownian Man, in this scenario the Lion has no 
winning strategy, while the Man is always sure to get away from the Lion. The interesting fact is that in 
this case the distance cannot grow exponentially fast, and in fact it has linear growth in time. To see this, 
from the previous inequality we obtain 

2arcsinh(e*'”'“^^*/^ sinh(/9(0)/2)) < p{t) < 2arccosh(e*'’^“^^*/^ cosh(p(0)/2)), 

which shows that p{t)/t converges to n — 1 as t —> oo, and therefore the distance function p{t) grows 
linearly in time for large times. 

Denoting by the model space of the n-dimensional manifold of constant curvature K, i.e. the 
Euclidean space for /F = 0, the sphere of radius 1 / ^/K for iF > 0, and the hyperbolic space in the case 
of FT < 0 (see Section 2 for the details), we obtain the following unified form of our results given in the 
abstract. This is the main result of this paper, and we include it here in a formal way. 


Theorem 1.1. For arbitrary distinct points x,y £ (n > 1) and a given non-negative function p : [0, oo) — )• 
[0, oo) with p{0) = d{x, y), there exists a co-adapted coupling of Brownian motions {X{t),Y (t)) on starting 
at {x,y) with deterministic distance function p{t) = d{X{t),Y{t)) if and only if p is continuous on [0,oo) and 
satisfies for almost every t >0the differential inequality 

(1.1) - - l)^/lf tan (=^) < m < - l)t/lf tan {^) + 

In particular, for dimension n = 1, the only co-adapted coupling {X{t),Y{t)) of Brownian motions on 
with deterministic distance is given by 


Y{t) 


'e*®X(t) if K>0 
< X{t) + e, ifK = 0 

ex{t), ifK<0 


for some 6 gR, with 6 > 0 if K < 0. 


Via a simple scaling argument, the proof of the theorem follows from the three cases discussed above: 
the case of the Euclidean space (Theorem 3.1), the case of the unit sphere (Theorem 4.1), and the case 
of the hyperbolic space (Theorem 5.1). 

The outline of the paper is the following. In Section 2 we introduce the basic notations and results 
needed in the sequel. For further use in the analysis of co-adapted couplings of spherical Brownian 
motions, in Lemma 2.1 we derive a characterization of all such couplings, similar to the one obtained 
in [11] or [15] in the case of Euclidean Brownian motions, and intimately related to Stroock's repre¬ 
sentation of spherical Brownian motion. Section 3 contains the analysis of Brownian couplings on M” 
(Theorem 3.1), and in Section 4 we present the analogous result for spherical Brownian motions (Theo¬ 
rem 4.1). In Section 5 we analyze the hyperbolic space case. The paper concludes with Section 6, which 
contains several remarks and corollaries of the main theorems, regarding the existence of fixed-distance, 
distance-increasing and distance-decreasing couplings, and an interpretation of the main results in terms 
of the stochastic Lion and Man problem. 

We point out that from the geometric point of view the construction of the couplings is an extrinsic one, 
relying on Stroock's representation of spherical Brownian motion in terms of a Brownian motion in the 
ambient Euclidean space and the realization of the hyperbolic space as the half space in Euclidean space. 
The advantage of this approach is that we can describe explicitly all co-adapted couplings, the downsize 
being that the construction applies only to model spaces. In another paper ([20]) we investigate, using an 
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intrinsic approach, partial extensions of the results obtained in the present paper to the case of smooth 
manifolds without boundary. 


2. Preliminaries 


We identify the vectors in with the corresponding n x 1 column matrices, and for a vector x G 
we denote by x' the transpose of x. The dot product of two vectors x,y £ M” will be written in terms of 
matrix multiplication as x ■ y = x'y, and we will denote by |x| = \/x'x the Euclidean length of the vector 
x G M”. We will use to denote the standard basis of W^. 

We denote the n-dimensional unit sphere in by 

§” = {x G : |x| = 1} 

and by d{x, y) the geodesic distance on it. The relationship between the geodesic distance on S” and the 
Euclidean distance is given by 

(2.1) d (x, y) = arcsin — {x'y)^ = 2 arcsin 

There are various ways of describing the spherical Brownian motion on S”, that is the Brownian mo¬ 
tion living on the surface of the sphere S” (see for example [6]). We will use Stroock's representation of 
spherical Brownian motion ([21]), as the solution X{t) of the Ito's stochastic differential equation 

(2.2) X{t)=X{0)+ [ {I - X{s)X{s)')dB{s)[ X{s)ds, t>0, 

Jo 2 Jq 

where B{t) is a (re -|- 1) -dimensional Brownian motion. The last term above may be thought as the pull 
needed in order to keep X{t) on the surface of the sphere In terms of Stratonovich integrals (see for 
instance [12, Chapter 2]), the above can be written equivalently as 

(2.3) X{t) = X{0)+ [ {I - X{s)X{sy) odB{s), t>0, 

Jo 

where "o" denotes the Stratonovich integration. Note that the operator I — X{t)X{t)' is simply the 
projection from onto the tangent space to the sphere at X{t). Equation (2.2) represents extrinsic 
formulation of the spherical Brownian motion on seen as a submanifold of 
The Eaplacian Ag on the sphere can be computed as 

{Asf)ix) = iAEf){x), xG§" 

where f{x) = / {x/\x\), and Ag is the Euclidean Eaplacian. Using (2.2) it is not hard to see that f{X{t)) — 
2 Jg AEf{X{s)) is a martingale for any smooth function / on S”, which proves that X{t) is indeed the 
spherical Brownian motion. 

We denote by S” = {x G : |x| = r} the re-dimensional sphere of radius r > 0. It is well known 
that the curvature of S” is constant equal to 1/r^, and for this reason S” is called the model space of constant 
curvature 1/r^. The Brownian motion X^{t) on S" is given by a simple rescaling of the Brownian motion 
X{t) on the unit sphere S”, and it can be described as X'^{t) = rX{tfr'^), with X{t) solving (2.2). 

There are several models for the manifold of constant curvature — 1, from which we choose here the 
half-space model. Our model for the n-dimensional hyperbolic space is thus given by 

H” = {(xi,...,x„) G M" : XI > 0}. 


(^) 


Eor a point x = (xi, X 2 ,..., x„) G we will often denote x = (0, X 2 ,..., x„). The metric on H” is 
given by 


2 dx\ -\- dx2 dx^ 


(2.4) ds^ = 

xf 

and the corresponding distance is given by 

(2.5) d(.,y) = arccosh ((l + " S''' 


2xiyi 


2xiyi 
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where arccosh denotes the inverse of the hyperbolic cosine function. 

The Laplace-Beltrami operator on (see for instance [10]) for details) is defined as 


( 2 . 6 ) 




i=2 


dxi' 


The Brownian motion on H” is the diffusion process on H” with generator given by A/2 above. A differ¬ 
ent way of constructing the process (see for example [18]) is as the solution of the stochastic differential 
equation 


(2.7) 


_2 

dX{t) = Xi{t)dB{t) --— Xi{t)eidt, 


where B{t) is a n-dimensional Euclidean Brownian motion. It is straightforward to check that the process 
X (t) defined in this way is indeed the Brownian motion on by simply using the Ito's formula and a 
straightforward calculation. 

The space defined above is the space of constant curvature —1. In a similar way we can define the 
model space of constant curvature — 1/r^ (r > 0) by e- = {x G : xi > 0}, with metric given by r'^ds'^, 
where ds‘^ is given by (2.4). The Brownian motion on this space is defined by X^{t) = rX{t/r‘^), where 
X{t) is the solution of (2.7). 

In order to give a unified statement of our results, we define the model space of dimension n > 1 and 
constant curvature A G M by 

if A > 0 




= 




P 


if A = 0 . 
if A < 0 

Recall that in general by a coupling we understand a pair of processes {X{t),Y{t)) defined on the same 
probability space, which are separately Markov, that is 

P {X{s + f) G A| A(s) = z, X{u) :0<u<s) = P^ {X{t) G A) 

P {Y{s + f) G A|y(s) = z,Yu : 0 < u < s) = P^ {Y{t) G A) 

for any s,t>0 and any measurable set A in the state space of the processes. 

The notion of Markovian coupling as used in [4] requires that in addition to the above, the joint process 
{X{t),Y(t)) is Markov and 

P {X{s + f) G A| A(s) = z, X{u),Yu : 0 < u < s) = P^ {X{t) G A) 

P (y(s + f) G A|y(s) = z, X{u),Yu :0<u<s) = P^ {Y{t) G A) 


( 2 . 8 ) 


for any s,t > 0 and any measurable set A in the state space of the processes. Markovian couplings 
are easily obtained for instance in the case when the process {X{t), Y (i)) is actually a diffusion. This 
would be the ideal case, but one can still construct a Markovian coupling by patching together diffusion 
processes for a conveniently chosen set of stopping times. 

The notion of co-adapted coupling (introduced by Kendall, [15]) is the same as the above but without the 
Markov property of {X{t),Y(t)). By a result on co-adapted couplings (Lemma 6 in [15]), a co-adapted 
coupling {X (t), Y (i)) of Brownian motions in can be represented as 

J{s)dX{s)+ [ K{s)dC{s), t>0, 

Jo 

where C is a n-dimensional Brownian motion independent of X (possibly on a larger filtration), and 
J,Ke M„xn are matrix-valued predictable random processes, satisfying the identity 

(2.10) J{t)J{ty + K{t)K{t)' = 1, t>0, 

with I denoting the n x n identity matrix. 

The formulae (2.9) and (2.10) provide an explicit representation of co-adapted couplings of n-dimensional 
Brownian motions, which will be used in order to characterize the couplings of deterministic distance in 
Euclidean spaces. In order to derive the equivalent characterization for couplings of spherical Brownian 


(2.9) 


Y{t) = y(o) + 
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motions, we need an explicit representation of co-adapted couplings on The following result gives 
such a characterization, intimately related to Stroock's representation (2.2) of the spherical Brownian 
motion. 


Lemma 2.1. Let {X{t),Y{t)) be a coupling of Brownian motions on S"- (n > 1). The coupling {X{t),Y{t)) is 
co-adapted if and only if there exists a co-adapted coupling W{t)) of [n -\-l)-dimensional Brownian motion, 

such that 


( 2 . 11 ) 


X(t)=x(0)+ [ {I - X{s)X{sy)dB{s) [ X{s)ds 
Jo 2 Jq 

y(f) = y(o)+ [ (/- y(s)y(s)')diy(s) - [ Y{s)ds 

Jo 2 Jo 


t > 0. 


Proof. Although the above is a statement about Brownian motions which is defined extrmsically in the 
language of differential geometry, in the proof we are going to use the intrinsic construction of the Brow¬ 
nian motion in terms of orthonormal frame bundle, as it is discussed for example in [12] and [22]. 

The converse implication is easier to prove, and we begin with it. Let X = (JJ)^^^ be the filtration 
generated by B and W. The hypothesis that the coupling {B{t), W{t)) is co-adapted shows that both 
B{t) and W{t) are Brownian motions with respect to the filtration F. Being solutions of (2.11), X{t) and 
y (f) are also adapted with respect to the filtration F. The conditions in (2.8) follows now easily from the 
fact that both X{t) and Y{t) are Markov processes with respect to the filtration F, hence they are also 
Markov processes with respect to the smaller filtration generated by (A, y). 

To prove the direct implication, assume that {X(t) ,Y(t)) is a co-adapted coupling of Brownian mo¬ 
tions on and let F = (JT)^>q be the filtration generated by {X, y). Condition (2.8) shows that both 
X and y are Markov processes with respect to the filtration F, and since they are also diffusions on S"^, 
with one half the Laplacian on as generator, it follows that X and Y are semimartingales with respect 
to the filtration F. To see this, considering a smooth and compactly supported function / : —)• M and 
denoting by Pt the semigroup generated by ^As on S”, we have 


Fs ] = iPt-sDiXis)) --J^ EiAsfiXiu))\Fs)du 


E(^f{X{t))-^l^ Asf{X{u))du 
= {Pt-sf){X{s)) Asf{X{u))du - ^ Pu-sAsf{X{u))du 

= {Pt-sf){X{s)) Asf{X{u))du - ^Pu-sf{X{u))du 

= {Pt.sf){X{s)) Asf{X{u))du - Pt-sfiXis)) + f{Xis)) 

= f{X{s))-\l\sf{X{u))du 


for all t > s > 0, which shows that f{X{t)) — ^ fQAsf{X{u))du is a martingale with respect to the 
filtration F. In turn this implies that A is a semimartingale with respect to the filtration F , with a similar 
proof for y. 

The second step in proving the existence of the Brownian motions B{T) and W{t) satisfying (2.11) is 
to use orthonormal frame bundles in order to construct their corresponding anti-development motions 
B{t) and W{t). The detailed construction is clearly presented in [12, Section 2.3], and we are not going 
to insist on it here. Note that by construction, B and W are JJ-adapted Brownian motions. 

Given the starting points x = A (0) and y = y (0) of the coupling, we can simply identify the tangent 
spaces TajS” and at these points with M”, and think of the Brownian motions B and W as moving in 
these spaces. To complete the construction, we consider 1-dimensional Brownian motions B{t) and W(t), 
independent of each other and also independent of Ft, which take values on the line determined by A (t) 
and Y{t). Let at be the parallel transport along the Brownian motion A from x to X{t), and /I* be the 
parallel transport from y to y (t) along the Brownian motion Y. It is easy to extend the operators at and 
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Pt to orthogonal matrices on by simply defining them to send the vector x into X{t), respectively 
y into Y(t) (note that both X(t) and Y (t) are in the space S"). 

Consider the decompositions x Mx, = Tj^S” x My (with Mx = {ax : a G M}), and 

define 

dB{t) = atd{B{t), B{t)) and dW{t) = Wt). 

Setting Qt to be be the ci-algebra generated by Bt and {Bg, : 0 < s < f}, it is a simple matter to 
check that B{t) and W{t) are ^j-adapted (n + 1)-dimensional Brownian motions. 

By construction (see [12, Theorem 2.3.4 and Lemma 2.3.2]) we have 

dX(t) = (I- X(t)X(ty) O (atdB(t)) 

and since {I — X{t)X{ty) {atdB{t)) = {I — X{t)X{ty) {atdB{t)), we obtain 

dX{t) = {I- X{t)X{ty) o dB{t), 

and similarly 

dY{t) = {I- Y{t)Y{ty) o dW{t). 

The fact that {B{t), W(t)) form a co-adapted coupling follows now easily, concluding the proof. □ 

We close this section with the analogous result given in the previous lemma, for the case of hyperbolic 
Brownian motion. 


Lemma 2.2. Let {X{t),Y{t)) be a coupling of Brownian motions on the hyperbolic space (n > 1). The 
coupling {X {t),Y (t)) is co-adapted if and only if there exists a co-adapted coupling {B{t), W (t)) of n-dimensional 
Brownian motions, such that 


( 2 . 12 ) 


ft ^ _ 2 /■* 

X{t)=X{0)+ / Xi{s)dB{s) -— / Xi{s)eids 


10 


/o 


ft „ _ 2 r* 

Y{t) = Y{Q)Y / Yi{s)dW{s) -— / Yi{s)eids 


t > 0. 


/o 
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Proof. The direct proof is rather simple. Let Xt be the sigma algebra generated by the Brownian motions 
B{s) and W (s) for s G [0, t]. Since the processes X{t) and Y(t) are defined by the stochastic differential 
equations (2.12), we see that the sigma algebra generated by Xg and Yg, s G [0, t] is included in Xt- Thus 
X (t) and Y (t) are Markov processes with respect to the filtration Xt, which proves the result. 

For the converse implication, consider the filtration X = {Xt)^^Q generated by {X, Y), and note that 
both X and Y are Markov processes with respect to the filtration X (by the co-adapted hypothesis), and 
therefore they are semimartingales with respect to X (the proof is similar to that in the previous lemma). 

Defining 

" L and W(t) = ^ yf-dW) + 


it is not difficult to check that B and W are Brownian motions with respect to the filtration X. For 
instance, we can use Levy's criterion and show that B{t) is an JT-nrartingale with quadratic variation 
given by tl, where I is the n x n identity matrix. The martingale characterization (2.7) of the hyperbolic 
Brownian motion X shows that the 2Li(t) — Xi{s)ds and X 2 {t), ..., Xn{t) are JT-nrartingales, from 

which we deduce that B{t) is also a JT-niartingale. Using again (2.7) and the definition of B above, we 
conclude that 

ft I 1 

{Bi,Bj)t = I 2 / ..d{Xi,Xj)g = I 2 / ^dijXi[s)ds = 6ijt, 1 ^ i, j ^ n. 

Jo ^iw) Jo 

A similar proof applies to W, and the claim follows. □ 
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3. Brownian couplings of deterministic distance in 

The main result of this section is the following characterization of all co-adapted couplings {X (t) ,Y (t) ) 
of n-dimensional Brownian motions, for which the distance |3r(t) — is a deterministic function of 
t > 0 . 

Theorem 3.1. For any distinct points x,y ^ R^ (n > 1) and an arbitrary non-negative function p : [0, oo) —)• 
[0, oo) with p{0) = \x — y\, there exists a co-adapted coupling of Brownian motions {X{t), Y (t)) in M” starting 
at {x,y) with deterministic distance function |X(t) — y(t)| = p{t) if and only if p is continuous on [0,oo) and 
satisfies almost everywhere the differential inequality 

(3.1) 0 < p'{t) < t > 0. 

p{t) 

In particular, the only co-adapted coupling of Brownian motions in M"- with (deterministic) non-increasing 
distance function is the translation coupling, and the only co-adapted coupling of Brownian motions in R with 
deterministic distance function is the translation coupling. 

Proof. If {X, Y) is an arbitrary co-adapted coupling of n-dimensional Brownian motions, Y can be rep¬ 
resented as in (2.9), where C is a n-dimensional Brownian motion independent of X, and the matrices 
J, K satisfy (2.10). 

Setting Z(t) = X(t) — Y(t) and using Ito's formula we obtain 

n 

d\Z{t)f =2Z{t)'dZ{t) + Y,d{Z^)t 

i=l 

n 

=2Z{ty (I - J{t)) dX{t) - 2Z{t)'K{t)dC{t) + ^ d{Z^)t. 

i=l 

Using the independence of X and C, and the relation (2.10) we obtain 

n 

d{Zyt = tr ((/ - J{t)) (I - J{t))' + K{t)K{ty) df = 2 (n - tr ( J(i))) dt. 

i=l 

From the last two equations we arrive at 

d\Z{t)f = 2Z{ty (I - J{t)) dX{t) - 2Z(ty K(t)dC(t) + 2 (n - tr (J(i))) dt. 

To prove the first claim of the theorem, note that if {X, U) is a co-adapted coupling of Brownian mo¬ 
tions with deterministic distance function p{t) = \Z{t)\, the martingale part of \Z{t)f must be identically 
zero, and therefore (using the independence of the Brownian motions X and C) we must have 

(3.2) J{tyz{t) = Z(t), K{tyz{t) = 0, and p^{t) = p‘^{0) + [ 2{n - tr {J{s)))ds 

Jo 

for all t > 0. 

From (2.10) it follows that 

0 < x' J{ty J{t)x < x' J(t)' J{t)x -|- x'K(tyK{t)x = x'x, X G M"", 

which shows that \ J{t)x\ < \x\, and in particular |e'J(f)ei| < 1 for z = 1,..., n, and thus — n < tr( J(t)) < 
n. Combining the above with the last equation in (3.2) we conclude that p is a non-decreasing function 
which satisfies 

(3.3) 0 < p{t)p'{t) = n — tr(J(t)) 

for almost every t > 0. In particular, this shows that p{t) > p(0) > 0, so Z{t) f 0 for any t > 0, and 
therefore we can find an orthonormal basis {^i(t),..., ^ri(i)} of such that ^i(f) = Using the first 
equation in (3.2), we obtain 

n n 

n-tr(J(t)) =tr(/-J(t)) = Y^i{t) {I - J (t)) Cft) = ^ (t)) Ci{t) ■ 

i=l i=2 
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Since = 1, we obtain J< 

|^j(t)| = 1. Combining with the last equation this shows that n— tr(J(t)) < 2(n — 1), which together with 
(3.3) concludes the proof of the direct implication. 

In order to prove the converse implication, we have to show that if the distance function p satisfies the 
equation (3.1) for almost every f > 0, we can find the matrices J{t) and K{t) satisfying (2.10) and (3.2). 

The hypothesis (3.1) shows that \Z{t)\ = p{t) is a non-decreasing function of f > 0, and in particular 
p{t) > p{0) > 0, which shows that Z{t) is a non-zero vector for any i > 0. Consider an orthonormal basis 
of M"" such that ^i(f) = and let J(t)'be the matrix with eigenvalues 1, X{t ),..., X{t) 
and corresponding eigenvectors ■ ■ ■, where X{t) G M is yet to be determined. 

With the above choice for J{t), we have J{t)'Z{t) = \Z{t)\J{t)'^i{t) = \Z{t)\^i{t) = Z{t), so the first 
equation in (3.2) is satisfied. We also have tr(J(t)) = 1 + {n — l)A(i), so the last equation in (3.2) reduces 
to p^{t) = /O^(0) -|- /g 2(n — 1) (1 — X{t)) dt, which is satisfied for the choice X{t) = 1 — ^^p(i)/0^(i) (and 0 
for definiteness if p'{t) is undefined, or if it does not satisfy (3.1)). Note that we have to choose X{t) only 
when n > 2, and in this case from the hypothesis (3.1) we have X{t) G [—1,1]. 

To complete the proof, we need to show that we can also choose the matrix K (t) such that the the sec¬ 
ond relation in (3.2) and the relation (2.10) are satisfied. Considering K(t) = PiDi ^P/, where Pt is the ma¬ 
trix with column vectors ■ ■ ■, and Di^t, ^ 2,4 are the diagonal matrices with diagonal en- 

triesO, y^l — A2(t),..., y^l — A2(t),respectively 1, A(f),..., A(f),wehavehr(f)'Z(t) = \Z{t)\PtD']^ ^Pl^l{t) 
0 and 

J{t)J{t)' + K{t)K{ty = Pt {Dl, + Dl,) Pi = PtPl = I, 

concluding the proof of the first part of the theorem. 

The second part follows now from the characterization given in the first part: if the distance function 
p in non-increasing, from (3.1) it follows that p'{t) = 0 for almost every t > 0, so tr(J(t)) = n by (3.2). 
The inequality (3.3) shows that the equality tr(J(i)) = n can hold iff J{t) is the identity matrix, so K{t) 
is the null matrix by (2.10), and therefore by (2.9) we obtain Y{t) = Y (0) — X(0) -|- X (t) for almost (hence 
all) t > 0, i.e. {X, y) is a translation coupling. 

Finally, the characterization given in the first part also shows that for n = 1 we have p' = 0 a.e., so 
p{t) = p(0) for all t > 0, and therefore {X, y) is a translation coupling, concluding the proof. □ 

4. Couplings of deterministic distance on S” 

In this section we study the co-adapted couplings of Brownian motions on the unit sphere 
The main result is the following. 

Theorem 4.1. For any points x,y € (n > 1) with x / ±i/ and an arbitrary non-negative function p : 
[0, cx)) —> [Ojoo) with p{Q) = d{x,y), there exists a co-adapted coupling of Brownian motions {X,Y) on S” 
starting at {x,y) with deterministic distance function d{X{t),Y{t)) = p{t) if and only if p is continuous on 
[0, CX)) and satisfies almost everywhere the differential inequality 

(4.1) 1) ton (^) </{()<-(«-!) tan t>0. 

In particular, the only co-adapted coupling of Brownian motions on with deterministic distance function is 
the rotation coupling, i.e. Y(t) = P^X{t)for some 0 G M and all t >0. 

Proof. To simplify the notation, we will prove an equivalent statement involving the function 

ri(t) = X'{t)Y (t) = cos {p{t)) 

Note that p is a deterministic function iff p is so, and the inequality (4.1) is equivalent to 

(4.2) — (n — 1)(?7 -|- 1) < p'{t) < —(n — l)(r/ — 1) 

Before proceeding to the proof, we will first argue that without loss of generality we may assume that 
X{t) f ±y(i) for all t > 0. Assume that the theorem has been proved for all set of times in t G [0, T), 
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where T = inf{t > 0 : r/(t) = 1 or rj{t) = —1}. The double inequality in (4.2) can be written in the 
equivalent form 

(4.3) — (n — l)rj{t) — (n — 1) < r/(t)' < —(n — l)'r]{t) + (n — 1), for a.e. t G [0, T), 

which combined with GronwalTs lemma (for which the almost ever}rwhere differentiability of rj suffices) 
gives 

-1 + (r?( 0 ) + l)e-("-^)* < r,{t) < 1 + ( 17 ( 0 ) - l)e-("-i)*, t G [ 0 ,T). 

If T were finite, by the continuity of r/ we would obtain 

{-1,1} 3 r 7 (r) C [-1 + ( 77 ( 0 ) + l)e-(-i)^, 1 + ( 77 ( 0 ) - l)e-(-i)^] C (-1,1), 

a contradiction. We conclude that without loss of generality we may assume that X (t) / ±.Y {t), for all 
i > 0, or equivalent that 77 (t) G (—1,1) for all t > 0. 

Assume now that (X, Y) is an arbitrary co-adapted coupling of Brownian motions on the sphere 
By Lemma 2.1 there exists a co-adapted coupling {B, W) of {n + l)-dimensional Brownian motions such 
that 

X{t) = A(0) + fl {I - X{s)X{s)') dB{s) - 5 X(s)ds 
Y(t) = Y(0} + (I - Y(s)Y(sy) dW(s) - ^ /o* Y(s)ds, 

and by Lemma 6 in [15] we also have 

J{s)dB{s) + f K{s)dC{s), 

Jo 

where C is a ( 77 - 1 - l)-dimensional Brownian motion independent of B, and J{t), K{t) are (n-|-1) x (n -|-1) 
matrices satisfying 

(4.5) + K{t)K{t)' = 1, t>0. 

Let U{t) = I — X{t)X{t)' and V{t) = I — Y{t)Y{ty (note that U{t) and V{t) are symmetric matrices, 
with U (f)^ = U (f) and V (i)^ = V(t)). We have 




n+1 

dT]{t) = X{tydY{t) + Y{tydX{t) + d{Xk, Yk)t 

k=l 

= {X{tyv{t)j{t) +Y{tyU{t))dB{t) + X'{t)V{t)K{t)dC{t) 

— nX{tyY{t)dt -|- ti{U{t)J{tyV{t))dt 

= {X{tyv{t)j{t) +Y{tyU{t))dB{t) + xyt)V{t)K{t)dC{t) 

— nr]{t)dt -|- ti{U{t)J{tyV{t))dt. 


Enforcing that 77 is a deterministic function requires the cancelation of the martingale part of it. Since the 
Brownian motions B and C are independent, this is equivalent to 

J{tyv{t)x{t) = -U{t)Y{t) and K{tyV{t)X{t) = 0. 

Using the definition of 77(f), U(t), and V (f), the above is also equivalent to 

(4.7) J{ty{X{t) - vit)Y{t)) = v{t)X{t) - Y{t) and K{ty{X{t) - v{t)Y{t)) = 0. 

To compute the trace ir{U {t)J{tyV (f)) in (4.6) we will use the following representation of the trace of 
a (77 -6 1 ) X (77 -6 1 ) matrix A: tr(A) = Ya=i where Afi = aijfj and /i,..., fn+i is an arbitrary 

basis of Since X{t) / ±Y{t), the subspace ff(f) spanned by X{t) and Y{t) is 2-dimensional, and 

let V(f) C be the (77 — 1)-dimensional subspace orthogonal to it. Using (4.7), we obtain 


u{t)j{tyv{t)x{t) = v{t)x{t) - Y{t) 
u{t)j{tyv{t)Y{t) = 0 . 
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If n = 1, X{t),Y(t) is a basis of U{t) = and from the above representation of the trace we obtain 
tv{U{t)J{tyV{t)) = which together with (4.6) proves the claim (4.2) in this case. If n > 2, consid¬ 
ering an arbitrary orthonormal basis /i(f),..., of Vt and using the above representation of the 

trace, we obtain 


n—1 

tr{u{t)j{tyv{t)) = vit) + Y,MtyU{t)JityV{t)Mt), 

i=i 

and combining with (4.6), we obtain 

n—1 

(4.8) d{r]{t)) = -{n - l)r]{t)dt + '^fj{tyU{t)J{tyV{t)fj{t)dt. 


The above proves the claim (4.2), if we show that each term in the sum above lies between —1 and 1. To 
show this, note that the condition (4.5) implies that | J(t)'C| < |^| for any vector ^ G so the operator 
norm of J(t)' is at most 1. It follows that the operator norm of J{t) is also at most 1, and since U{t) and 
V{t) are projection operators, it is trivial to see that |/,(f)'17(f)J(f)'l/(f)/j(f)| < l,foranyj = l,...,n — 1, 
concluding the proof of the direct implication of the theorem. 

To prove the converse, we have to show that if rj satisfies (4.2), we can find (n -|- 1) x (n -|- 1) matrices 
J{t),K{t) which satisfy (4.5), (4.7) and (4.8). To this aim consider the matrix J{t) defined by 

'j{tyx{t) = -Y{t) 

< j(t)'y(f) = -x{t) 

J{tyV = 7(t)y, for any V G V(f) 

where 

= 1 _ + (n - l)v{t) 

n — 1 ’ 

which from the condition (4.2) is guaranteed to belong to [—1,1]. Consider also the matrix K{t) defined 
by 

' K{tyx{t) = 0 
< K{tyY{t) = 0 

K{tyV = \/l — for any V G V{t) 

and note that with the above choices we have J{t) = J{ty, K{t) = K{ty, J{t)J{ty + K{t)K{ty = I. In 
addition it is easy to check (4.7) and (4.8), hence concluding the proof. □ 


5. Couplings of deterministic distance on the hyperbolic space 


Theorem 5.1. For any distinct -points x,y e 'EY (n > 1) and an arbitrary non-negative function p : [0, oo) — )• 
[0, oo) with p{0) = d{x,y), there exists a co-adapted coupling of Brownian motions {X,Y) on starting at 
{x,y) with deterministic distance function p{t) = d{X{t),Y{t)) if and only if p is continuous on [0,oo) and 
satisfies almost everywhere the differential inequality 

(5.1) t > 0. 

In particular, in dimension n = 1, the only co-adapted Brownian coupling {X, Y) of deterministic distance on 
is given by Y(t) = OX(t), t > 0,/or some 9 > 0. 


Proof. Before proceeding with the proof, note that by Lemma 2.2, the Brownian coupling {X, Y) on H” is 
co-adapted iff there exists a co-adapted coupling {B, 'W) of Euclidean Brownian motions such that 

X(t) = X(0)+ / Xi{s)dBs -— / Xi{s)eids 


10 


/o 


/■* n — 2 /■* 

Y{t) = y(0) + / Yi{s)dWs - / Yi{s)eids 

Jo ^ Jo 


t > 0, 
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and dW{t) = J{t)dB{t)+K{t)dC{t), with B, C independent n-dimensional Euclidean Brownian motions 
and J{t),K{t) n x n matrices satisfying + K{t)K{t)' = I. Denoting Z{t) = X{t) — Y{t), the 

statement of the theorem can be expressed equivalently in terms of the function 


(5.2) 


d{t) 


\z{t)? 

2Xi(t)yi(t)' 


The equation (2.5) shows that cosh(d(X(t), D(t))) = r?(t) + 1, and the inequality (5.1) is therefore equiva¬ 
lent to 


(5.3) f?(i)(ra “ 1) < d'{t) < “ 1) + 2(n — 1), t > 0. 

In order to prove the direct implication, we assume that {X, y) is a co-adapted Brownian coupling on 
with deterministic distance p and show that (5.3) holds. We first note, just as in the spherical case, 
that without loss of generality we may assume that X{t) / Y{t) for all t > 0: if (5.3) is satisfied for 
0 < f < T = inffs > 0 : 7 y(s) = 0}, then 

7?(0)e("-i)* < 77(f), fG[ 0 ,r), 

and thus rj does not vanish on [0, T). If T were finite, by the continuity of p we would obtain 77 (T) = 0, a 
contradiction. Consequently we may assume that X{t) i^Y (t) for all t >0. 

Using the preliminary remarks, we obtain 

Th - 2 

dZ{t) = iXi{t)I - Yi{t)J{t))dB{t) - Yi{t)K{t)dC{t) - —Zi{t)eidt, 

and since p{t) is a function of bounded variation, from (5.2) we obtain 

(5.4) d\Z{t)\‘^ = 2Xi{t)Yi{t)p'{t)dt + 2p{t)d{Xi{t)Yi{t)). 

Since 

d{Xi{t)Yiit)) = Yi{t)e[dX{t) + Xi{t)e[dYit) + d{Xi,Yi)t 

= Xi{t)Yi{t)e[ (/ + J{t)) dB{t) + Xi{t)Yi{t)e[K{t)dC{t) 

+ Xi{t)Yi{t) {e[J{t)ei - {n- 2)) dt, 

combining with the above we conclude 

d\Z{t)\‘^ = 2p{t)Xi{t)Yi{t)e[ (/ + J{t)) dB{t) + 2p{t)Xi{t)Yi{t)e[K{t)dC{t) 

+ 2Xi{t)Yi{t) (r]{t) (e'iJ(f)ei - (n - 2)) +p'{t)) dt, 

Alternately, we can determine the semimartingale decomposition of |Z(f)p as follows: 

n 

d|Z(f)|2 = 2Z{t)'dZ{t) + J2d{Zi)t 

i=l 

n 

= 2M{t)dB{t) + 2N{t)dC{t) - (n - 2)Z^(t)dt + ^ d{Zi}t, 

i=l 

where 

M{t) = X,{t)Z{ty - Yi{t)Z{tyj{t), N{t) = -Yi{t)Z{tyK{t), 


i=l 


tr {{Xi{t)I - Yi{t)J{t)) {Xi{t)I - Yi{t)J{t)y + Yy{t)K{t)K{ty) dt 
[nXf{t) + nYy{t) — 2Xi{t)Yi{t)ti{J{t))) dt. 


and 
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Comparing the two semimartingale representations of |Zp above, and using the independence of the 
Brownian motions B and C, we conclude 


' M(t) = r]{t)Xiit)Yi{t)e[iI + J{t)) 

< N{t) = 7j{t)Xi{t)Yi{t)e[K{t) 

^ Xf{t) + Y^{t) + Xi{t)Yi{t) ((1 + 77(f)) (n - 2 - e'iJ(f)ei) - -tr(J(f))) = 0. 

where tr(J(f)) = Using the definitions of M{t), N{t), and 77 (also recall the notation 

X = { 0 ,X 2 , ■ ■ ■, Xn)), the above conditions are equivalent to 

' J{ty [{imy + xUt) - Y^it)) ei + 2Y,{t)Z{t)) 

= (^{-\Z{t)y + Xf{t)-Yy{t))e^ + 2X,{t)Z{t)) 


and 

(5.6) 


K{ty 

77 '(f) = 



+ xf{t)-Yyit)) 


ei + 2Yi 


(1 + 77 (f)) (n - 2 - e[J(t)ei) - tr(J(f)) + 


= 0 


xHt) + Yy(t) 

Xi(t)Yi(t) 


Consider first the case when Z(t) = 0. Since Z(t) / 0 we must have Xi(f) / Ui(f), and from the first 
equation in (5.5) we obtain J(f)'ei = ei. Using this, 1 + 7](t) = 2 x)(t)y(( 7 ) ' gives 

(5.7) rf'(t) = r]{t){n - 1) + (tt - 1) - tr(J(t)). 

Since J{t)J{ty + K{t)K{ty = I, we obtain —1 < e[J{tyei < 1, so —{n — 1) < tr(J(f)) < n — 1, and thus 
we arrive at (5.3). 

For the case of Z{t) 7 ^ 0 we will use the following result. 

Lemma 5.2. Assume A : M"' —)• M”' is a linear map such that \ Ax\ < |a:|/or any x G M”. Let ^ 1,^2 £ be 
orthonormal vectors so that 


( 5 . 8 ) , 4 ( 777,^1 + 1 ^ 2 ) = P^i + q^2 

for some m,l,p,q eR with mf + P = f 0. The following hold true: 

( 1 ) iIn-A'A)y^m^i+l^2) = 0. 

(2) We have for any r, s G M, 

rmp + slq — \rlq + smp\ rmp + slq + \rlq + smp\ 

- < r^iALi + 8^2 A^2 < - 


777 ^ + P 

(3) Moreover, for each r, s, (/> G M satisfying 

rmp + slq — \rlq + S 777 ,p| 


(5.9) 


777,' 


+ P 


< 6< 


m? + P 
rmp + slq + \rlq + smp\ 


m- 


+ P 


we have + s^ 2 A<pC 2 = 4>f where A^f, is the map defined by 

mp + dlq mq — dip 

"‘♦fi = 

/ , „ Ip — dmq _ Iq + dmp _ 
M 2 = + I 2 , ,9 6 




+p 


m‘ 


+ p 


^A,pi = 0 , for i orthogonal to ^ 1 , 6 , 


with d = 

rlq+smp 


Proof. (1) Note that if U is a symmetric, non-negative definite matrix for which = 0, then \ D^Ppf = 
= 0, and therefore = 0. Using this with D = I — A' A (which is non-negative definite) 

and ^ = 777,^1 -|- 1^2, one can check by direct calculation that = mp + P — p^ — q^ = Q from 
which the first claim follows. 
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(2) It is clear that we have to focus on the action of the matrix A on the space generated by and 
^ 2 - For convenience, let V be the subspace orthogonal to and consider the restriction of A to 
the subspace ZY, more precisely, A^ = 1174^ for where n is the orthogonal projection onto 

Writing vectors in ZY in terms of the orthonormal basis , ^ 2 , we can simply assume that MA is 
replaced by M^, and the map A is replaced by A. It is clear that the relation (5.8) continues to hold 
with A replaced by A. For simplicity we will also assume that m? + P = = 1. Consider 


We obviously have 


U 


m —I 
I m 


U 



m 

I 


so AUei = Vei. Denoting B = V'AU = 


and V 


p -q 
q p 


and V 

a b 
c d 


■ 1 ■ 


p 

0 


_ Q _ 


, then 


Bei = ei 


which means that o = 1 and c = 0. On the other hand, we also have that \Bx\ < \x\ for any 
X G M?, which implies that for x = ^ , 

{t + bs)‘^ + dFs^ <P A 


for any real numbers t,s. Thus, we must have 6 = 0 and — 1 < / < 1. Therefore, 


A = VBU' 


mp + dlq pi — dmq 
qm — dpi Iq + dmp 


Thus, 

re'iAei + 562^62 = rmp + slq + d{rlq + smp). 

Since — 1 < cZ < 1, the rest follows easily 
(3) This follows from a direct calculation and is left to the reader. 


□ 


Using the above lemma with A = J(t)', .^1 = ei, ^2 = Zl\Z\,m = \Zt\‘^ + X^t) — I = 2Yi{t)\Zt\, 

p = + Xf{t) - Y^{t), q = 2Xi{t)\Zt\, r = 1 + r]{t), and s = 1, we have 

rmp + slq X‘}{t)+Y?(t) ,, , rlq + smp 

-^- dV’) - 1 and -^= 1, 

and therefore we obtain 


where 


Xi{t)Yi{t) ^ m2 + Z2 

a - 1 < (1 + rj{t))e[J {ty ei + ^ 2 <F(Z )''^2 < a + 1, 

xf{t) + Yy{t) 


a = 


-r]{t) - 1. 


Xi{t)Yi{t) 

Completing {ei , ^ 2 } to an orthonormal basis {ei, .^ 2 , • • • , ^n} of M", and recalling that — 1 < Pidity^i < 1 
for any unit vector ^j, we have 

a — n + l<(l + p{t))e[J {t)' ei + tr( J(Z)) < a + n — 1, 

and therefore using (5.6) we obtain 

r]{t){n — 1) + (n — 2) — tr(J) < rj'{t) < ri{t){n — 1) + n — tr(J), 


concluding the proof of the direct implication. 

Conversely, assuming (5.3) holds, we have to show that we can choose matrices J{t),K{t) such that 
J{t)J{ty + K{t)K{ty = I, (5.5) and (5.6) are satisfied. 

In the case n = 1, choosing J{t) = —1 and K{t) = 0 proves the claim. 
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Consider now the case n > 2. If Z{t) = 0, the equation (5.6) is equivalent to (5.7), and therefore the 
claim is satisfied for the choice 


J{t)ei = 




i = 1 
i / 1 


and K{t)ei = 


0 i = l 

yjl- -i‘^{t)ei i / 1 


where 

(5.10) 


7(*) = l + t((i)-^e[-l,ll. 
n — 1 


For the case Z{t) ^ 0, we will use Lemma 5.2 as follows. Consider = ei, ^2 = m = 

\Z{t)\^ + Xl{t)-Y^{t),l = 2Yi{t)\Z{t)\,p = -\Z{t)\^ + Xl {t) -Y^{t),q = 2Xi{t)\Z{t)\,r = l + r^{t), and 
s = 1 , for which 


rmp + slq _ X^(f) + Yi{t) 


m‘ 


+ P 


Xi{t)Yi{t) 


— r]{t) — 1 and 


rlq Y smp 


= 1 . 


ni' 


Using the notation of the lemma, define the n x n matrices J and K by 




7(t)C, i G v(f) 




e G v{t) 


where U{t) = spanjei, Z{t)}, V{t) C M” is the subspace orthogonal to it, 7 (f) is given by (5.10), and 
(whose value will be chosen below) satisfies 

Xi(l)Vi(l) Xi(t)Vi(l) 

Note that I — is a symmetric non-negative definite matrix and 7 (f) G [—1,1], so the matrix 

K{t) is well defined (for definiteness, if ri'{t) is not defined, or if it does not satisfy (5.3), we consider 
7(f) = 0). 

With the above choice for J and K, the condition J{t)J{ty -|- K{t)K{t)' = I is clearly satisfied for any 
4>{t) satisfying (5.11). The first equation in (5.5)can be written equivalently as 

J{t)' (m^i +1^2) = (m^i +1^2) =pii + qi2, 

and is satisfied by the definition (5.9) of the matrix for any satisfying (5.11). The above relation 
(together with | 7 (f) | < 1) also shows that we can apply the first part of Lemma 5.2 with A = A^(^p thus 
obtaining 

K'(t) (m^i -6 1 ^ 2 ) = {I - + ^^ 2 ) = 0 , 

so the second equation in (5.5) is also satisfied, for any any verifying (5.11). 

Finally, we will choose (/)(f) so that the equation (5.6) is also satisfied. The last part of the lemma shows 
that 


(1 -6 7(f))e'iJ(f)ei -6 = </>(*)> 


7(f\' 

and completing {ei, jhpA to an orthonormal basis of M", as in the previous part of the proof, we obtain 
(1 -|- p{t))e'iJ{t)ei + tr(J(f) = 4>{t) + {n — 2)'y{t). The condition (5.6) is thus equivalent to 

xfit) + Yyit) 


rj\t) = {n-2){l + 7 (f) - 7 (f)) - (/.(f) -h 

and recalling the choice (5.10) of 7 (f), we are led to the choice 


Xi(f)yi(f) 


zm - ^ xf{t) + Ynt) 

n-C Xi(f)yi(f) • 

For the above choice of (/>(f) the condition (5.11) is satisfied (it is just the hypothesis (5.3)), thus con¬ 


cluding the proof of the theorem. 


□ 
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6. Corollaries and Remarks 
Several remarks and corollaries are in order here. 

Remark 6.1. In Theorem 1.1 we obtained a characterization of all Brownian co-adapted couplings for 
which the distance function p is deterministic. If the distance function p is in addition a solution to an 
ordinary differential equation of the form 

p'it) = F{p{t)) 

where F : [0, oo) —M is a function such that 

— (n — l)y/Ktsin{y/Kp/2) < F{p) < —(n — 1)Vk tan(\/iT/o/2) + ■ , 

s\n.{y Kp) 

the construction given in the theorem also provides couplings {X, Y) which are diffusions on the corre¬ 
sponding spaces. The argument is simply the observation that the construction of the matrices J{t),K{t) 
in this theorem depend only on X {t) and Y {t) and not on the time t > 0. 

Considering the case of the constant distance function p in Theorem 4.1, we obtain the following 
important result. 

Corollary 6.2. For any points x,y € (re > 2) there exists a co-adapted fixed-distance coupling, that is there 
exists a coupling {X, Y) of Brownian motions on starting at {x, y)for which d{X{t),Y (t)) = d{x, y)for all 
t > 0. Moreover, the joint processes {X, Y) is a diffusion. 

Remark 6.3. After writing this article, the second author learned from a private discussion with Thierry 
Levy an alternate short proof of the corollary above in the case re = 2, which we briefly present below. 
Consider 0(3), the set of 3 x 3 orthogonal matrices, and for each x G consider the map : 0(3) —)> 
given by = 2\.x. Consider the standard left-right invariant metric on 0(3), and the Riemannian 

structure associated to it. Denoting by and the Laplacian on 0(3), respectively on S^, it can 

be shown that o vr^) = A^^(/), and as a consequence, if Z{t) is a Brownian motion on 0(3), 

then Z{t)x is a Brownian motion on §^. For arbitrary points G S^, choosing X{t) = Z{t)x and 
Y(t) = Z{t)y, provides a co-adapted fixed-distance coupling {X{t),Y (t)) of Brownian motions on 

Considering p{t) = e~^^^‘^p{0), with A; < re — 1 in Theorem 4.1 we obtain the following. 

Corollary 6.4. For any points x, y G (re > 1) with x / ±y and 0 < k < n — 1, there exists a co-adapted 
coupling {X{t), Y (t)) of Brownian motions on such that d{X{t),Y (t)) = e~^^/‘^d{x, y)for all t >0. 

For k < 0 and re > 1, there also exists such a coupling, but it only satisfies d{X{t),Y{t)) = e~^^^‘^d{x,y), 
0 <t < to, for some to < oo. 

Note that in the case A; < 0 one cannot hope that the equality d{X{t),Y(t)) = e“^*/^d(x, y) holds true 
for all times f > 0 , since the exponential term is unbounded and the sphere is compact. 

An interesting feature of the previous corollary is the fact that the maximum allowable value of k is 
re — 1. It turns out that this has to do with the fact that the curvature of the sphere S" is 1 (more precisely 
the lower bound on the Ricci curvature is re — 1, but the reasoning is not transparent from the extrinsic 
argument presented in the proof of the theorem. 

Remark 6.5. Aside from the particular case when the function p is constant, there are two other interesting 
particular cases of Theorem 4.1. One is the extreme case in which the left inequality in the hypothesis 

(4.1) of the theorem is attained, namely 

(6.1) p{t) = 2 arcsin sin(/ 9 ( 0 )/ 2 )^ , A > 0 . 

The corresponding coupling is a particular case of shy coupling, in which the two processes do not 
couple, but approach each other exponentially fast. 

The other case is the extreme case when the right side of the inequality (4.1) is attained, explicitly 

p{t) = 2 arccos cos(/ 9 ( 0 )/ 2 )^ , 


( 6 . 2 ) 


t > 0 . 



COUPLINGS OF BROWNIAN MOTIONS OF DETERMINISTIC DISTANCE 


17 


The corresponding coupling is again a shy coupling (the processes do not couple in finite time), but it is 
a repulsive coupling, in the sense that the distance between the two processes increases at an exponential 
rate to the maximum distance allowed on the sphere (the processes become antipodal in the limit). 

Notice that the latter coupling is related to the former coupling via the following simple observation. 
If {X, Y) is the coupling for which the distance function is given by (6.1), then (X(t), Y(t)) with Y(t) = 
—Y(t) is a coupling for which the distance function is given by (6.2). 


Remark 6.6. For the case of the hyperbolic space 'EY, we point out that for any function p : [0, oo) —> [0, oo) 
for which there exists a co-adapted coupling of Brownian motions on El” with d{X{t),Y (t)) = p{t), from 
(5.1) we get that 

2 arcsinh(e^”“^^*/^ sinh(/j(0)/2)) < p{t) < 2arccosh(e^”“^^*/^ cosh(/9(0)/2)), 


which in turn shows that 


lim 

t^OO 


pjt) 

t 


= n — 1. 


In particular, this shows that there cannot be an exponential growth at infinity for the distance function. 

Moreover, any function p{t) which can be realized as the distance between co-adapted Brownian mo¬ 
tions must be increasing. In particular, there are no fixed-distance couplings as in the case of the sphere 
or of the Euclidean space. 

It is also interesting to point out that for p{t) = e~^^p{Y), the inequality in (5.1) is satisfied only for 
/c < 0 and for small t. In turn, if we take p{Q) sufficiently small, we conclude that k < —{n — 1) /2. 


Though we studied the coupling of Brownian motions on the model spaces, we can extend part of 
the analysis to the case of constant curvature manifolds. To state this, let M be a manifold of con¬ 
stant curvature K and assume that i{M) is the injectivity radius of M. If p : [0,T) —^ [0,i(M)) is a 
continuous function such that (1.1) then we can find a coupling of Brownian motions {X, Y) such that 
p{t) = d{X{t),Y(t )). The argument is simply based on the fact that M is the quotient of by a discrete 
group of isometries. Since the Brownian motions on stay within i{M) distance, it turns out that their 

projections onto M are Brownian motions staying within the cut locus of each other. 

The initial motivation for writing the present article was to investigate the stochastic version of the 
Lion and Man problem presented in the introduction, and we conclude with an interpretation of our 
results in terms of it. Given a distance function p satisfying the hypotheses of Theorem 1.1, the Brownian 
Man can find a "strategy" to keep the Brownian Lion at distance p{t) at time t. In particular the Man can 
always find a strategy which keeps the Lion at fixed distance in the Euclidean space and on the sphere, 
which although it is trivial in the Euclidean case, it is completely non-trivial on the sphere. The Man can 
also find a strategy which increases his distance from the Lion in all model spaces. 

While the Euclidean and the hyperbolic cases must be disappointing for the Lion (no distance-decreasing 
coupling), in the case of the sphere the Lion can find a strategy which brings him exponentially close in 
time to the Man, which should be sufficient for some practical purposes. 

Acknowledgements. We want to thank Wilfrid Kendall and Krzysztof Burdzy for several interesting 
discussions on the existence of fixed-distance coupling on the sphere which took place in the summer 
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and complete this program on manifolds, which will appear in a forthcoming paper. The second author 
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